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It is shown that the spontaneous magnetization occurs due to the anomalous magnetic 
moments of quarks in the high-density quark matter under the tensor-type four-point interac¬ 
tion. The spin polarized condensate for each flavor of quark appears at high baryon density, 
which leads to the spontaneous magnetization due to the anomalous magnetic moments of 
quarks. The implications to the strong magnetic held in the compact stars is discussed. 

§1. Introduction 

One of the recent interests about the physics governed by quantum chromody¬ 
namics (QCD) may be to clarify the phase structure of QCD with respect to the 
temperature, baryon chemical potential or baryon number density, external mag¬ 
netic field and so on.® In the region with high temperature and low density, the 
lattice QCD simulation gives many insights about the QCD world. However, as is 
well known, in the region of large quark chemical potential, there are some prob¬ 
lems to be solved in order to calculate the physical quantities definitely. But, it 
is expected that there are still many interesting phenomena in the region of low 
temperature and large chemical potential. Especially, it is expected that there are 
many exotic phases such as the two-flavor color superconducting phase, the color- 
flavor locked phase,®’®’® the quarkyonic phase,® the phase with the inhomogeneous 
condensate and so forth. 

In the heavy-ion collision experiments such as relativistic heavy-ion collider 
(RHIC) experiments at Brookhaven National Laboratory, it is believed that the 
quark-gluon phase is realized apart from the hadronic phase. Thus, the quark mat¬ 
ter may be created and the extreme states of QCD with finite temperature and 
density may be realized in the heavy-ion collision experiments. On the other hand, 
the high-density hadronic phase or quark phase may be realized in the inner core of 
the compact star objects such as the neutron stars, magnetars and quark stars, if 
they exist. In the core of these compact stars, it is expected that the hadron or quark 
phase exists at low temperature and high baryon density. For example, the quark 
matter is considered to exist in the core of pulsars with two solar masses in compact 
stars, and in the case of the heavy ion collisions.® Therefore, the investigation of 
quark matter in the region of the low temperature and the large quark chemical 
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potential is one of very interesting and important subjects in order to understand 
the whole world governed by QCD. 

It is further known that neutron stars, especially the so-called magnetars,®® 
have a strong magnetic field. However, the origin of the strong magnetic field is not 
so clear. It has been pointed out that the strong magnetic field may be created if 
the quark liquid exists in the core of the compact stars.® Further, the possibility 
of the quark spin polarization in the high-density quark matter has been investi¬ 
gated when the pseudovector-type interaction between quarks exists.® Under the 
pseudovector-type interaction, it was shown that the spin polarized condensate ap¬ 
pears, which leads to a ferromagnetism in a quark matter.® However, the spin 
polarized condensate only appears in a narrow region of the quark chemical poten¬ 
tial. 

As for the spin polarized condensate, it has been shown that the tensor-type 
interaction between quarks leads to the spin polarized phase in high-density quark 
matter.®’®’® The present authors have shown the possibility of the quark spin 
polarized phase in the quark matter at high baryon density against the two-flavor 
color superconducting phas^® and the color-flavor locked phas^® under both the 
quark pairing interaction and tensor-type four-point interaction in the Nambu-Jona- 
Lasinio (NJL)-type model.®’® However, the magnetic feature is not considered in 
the tensor-type interaction until now . 

In this paper, the magnetic features are investigated under the existence of 
the spin polarized condensate in the NJL model with the tensor-type four-point 
interactions. It will be shown that the anomalous magnetic moments of quark^® 
play an essential role, the existence of which leads to the spontaneous magnetization 
of quark matter. Also, the implication to the magnetic field of compact stars such 
as neutron stars with quark matter, namely hybrid stars, will be discussed. 

This paper is organized as follows: In the next two sections, §2 and §3, the 
thermodynamic potential under the external magnetic field is given with and without 
the spin polarized condensate. In section §4, an approximate expression of the 
thermodynamic potential with small magnetic field is derived. In section §5, it 
is shown that the spontaneous magnetization does not appear in the case of no 
anomalous magnetic moments of quarks. In section §6, the anomalous magnetic 
moments of quarks are introduced within the mean field approximation. As a result, 
spontaneous magnetization occurs in quark matter in the region of high baryon 
density due to both the spontaneous spin polarization and the anomalous magnetic 
moments of quarks. In section §7, the implication to the hybrid compact stars is 
discussed briefly and it is shown that the strong magnetic held in the surface of 
the compact stars can be revealed. The last section is devoted to a summary and 
concluding remarks. 
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§2. Thermodynamic potential under external magnetic field 


We consider the two-flavor case. Let us start from the following Lagrangian 
density with chiral symmetry: 

G - 

To = ('07/^71^+ (0*757^7'''0)(0n57Ai7i/'0)) , (2-1) 

where r represents the flavor su( 2 )-generator. 

In previous study, we have found that a spin polarized condensate 
L 3 = may be realized at high baryon density, where Lio = —^ 7 ^ 7 ^ 

is the spin operator. Thus, the Lagrangian density under the mean field approxima¬ 
tion is obtained as 


Tmf = - — , 


( 2 . 2 ) 


T 3 = -G(0L'3r30) = Ftj , 


*12 

= -^7 7 = 


0-3 0 

0 as 


where tj = 1 for u-quark (/ = u) and tj = —1 for d-quark (/ = d), denote the 
eigenvalues of r 3 . Here, <73 is the third component of the Pauli spin matrices. 

Hereafter, let us consider the system under the external magnetic field B along 
the z-axis. The Lagrangian density is recast into 


£ = - Fs'ipFsTs'il) - — , 


where Dn represents the covariant derivative: 


Dn — da F iQA 




Aa={ 0 , 


By Bx 


0 = (0, -A) . 


(2.3) 


(2.4) 


Here, for up (down) quark, Q = 2e/3 (—e/3) where e is the elementary charge. 
Because we investigate the quark matter at finite density system, the Hamiltonian 
density with quark chemical potential y, is obtained as 

B. - yN = ipi-ij ■ (V - iQA) - ^7°)0 + FsipFsTsip + — 

= 0^(/i-/i)V’ +^ , (2.5) 

where Af represents the quark-number operator. Here, h can be expressed as 

h = a- {p- QA) + Fsts^Bs 

= azPz + oixPx + aypy + FsTsjSUs , ( 2 . 6 ) 

where, by using the Dirac representation of the Dirac gamma matrices. 


a* = 7°7* = 


0 cJi 

cj,' 0 


0 = 7 = 


1 0 
0 -1 


Pi = 


A 


Px=Px + , 


P ^ QB 

Py=Py- —X ■ 


(2.7) 


Hereafter, let us consider two cases, namely F > 0 and F = 0, respectively. 
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2.1. F > 0 case 

The Dirac equation is written by id'll) jdt = hij; = Eil), namely 

f FsTsaz - + PxC^x + Py(yy W A g ^2 8) 

V PzCTz + PxC^x + PyCTy -F^T^a^ - E J \ if J 

where 0 and ip are two-component spinors. Eliminating ip and noting \ Px , Py]~ 
iQB, we obtain the following equation: 



2F 



(2.9) 



for up quark 
for down quark 


where 1 is the identity matrix in the isospin space. Let us consider Q > 0 case, that 
is, the case of up quark. We introduce new operators instead of Px and Py as 


1 




( 2 . 10 ) 


where the commutation relation [ a , ] = 1 is satisfied. Then, Eq. (| 2.9I) is recast 

into 



( 2 . 11 ) 


Here, expressing the two-component spinor cp as 4> = the above equation 

can be expressed as 


Xo + To + (^1 + Yi)a^a <Pi - Ca)<P2 = 0 , 
-CacPi + [Xo - To + (^1 - Yi)a^a\ P)2 = 0 


( 2 . 12 ) 


and eliminating 02, the following equation can be derived from Ea. ()2.11l) : 

[Xo + To + (^1 + Yi)a^a] - [Xq - Tq + (^i - Yi){a^a - l)]-^C\^a\ 0i = 0 , 

(2.13) 
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where we define the following quantity for simplicity: 


Xo = -E + 


Epl 


+ 


E-F 


E^ - F'^ ^2 _ p2 


QB 


Xi = 


2EQB 
E"^ - F 2 


Id = 


2FQB 

^2 _ i ?2 


Yo = F- 
C = 


Fp 


E-F 

p-i- F‘^~ E"^ - F2 
2Fp^y/2QB 


E^ _ i ?2 


QB , 
(2.14) 


Here, pz and a^a should be replaced by their eigenvalues pz and v (= 0,1, 2, • • •). 
Further, since the case Q > 0 is treated, namely the case of up quark, so F 3 = F. 
Then from the coefficient of (pi being 0, we can get the eigenvalue of the Dirac 


^flavor 


as 


equation, E, which is expressed as j, , 

= {F ± i/2^Y + pI (for 1 / = 1, 2, • • •) 


gUP 

Pz,V 


(2.15) 


with <5 = 2e/3. It should be here noted that it is necessary to pay a special attention 
for the case 12 = 0 because a(p{v = 0) = a\v = Q) =0 satishes. For z/ = 0, from 
(pi = |z^ = 0), Ea. ()2.11l) is reduced to 


-^0 + ^0 — 0 , <?i>2 — 0 


(2.16) 


because acpi = a|i^ = 0) = 0. Form Xq + I 0 = 0) which leads to a quadratic equation 
with respect to E, we have 


up 

^Pz,U=0 




(2.17) 


for positive energy solution. Thus, for = 0, the solution only appears once, and 
the energy corresponding to z/ = 0 is not degenerate. 

In the same way, we can get the eigenvalue E in the case of down quark where 
Q < 0. Instead of () 2 . 10 p . we define new operators as 


1 


a = 


:(He - iPy) , 


a) = 


1 


]{Px + ‘^Py) ) 


(2.18) 


where [ a , ] = 1 is satisfied. Then, we obtain the equation instead of (12.1111 as 

E , , 4 . + ^ F 


-E + 


+ <^z i P — 


E2-F2 
F 


{pI — QB{aa^ + a^a) — 


E‘i - f 2 


{pI + QB{aa^ + a^a) — 


E^ - F 2 
E 


E‘2 - i ?2 


QB 
QB 


2F 


E^ - F 2 


Pz 


+ a^) + iay{a - a^)) 


= 0 . 


(2.19) 


Eliminating cpi which is the upper component of the two-component spinor cp, we 
obtain the following equation: 

[Xo + To - (^1 + Yi)a^a] - [Xq - - (Xi - Yi)(ata - l)]-^C\^a] (P2 = 0 , 

( 2 . 20 ) 
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where 
Xq = -E + 
Xi = 


Epl 


E^ - F 2 

2EQB 


QB , 


E + F 
£■2 - F 2 
2FQB 


E'i -F"^ ’ 


5^1 = - 


E^-F"^ ’ 


Yo = -F + 
C = - 


Fpl 


+ ^r^QB , 


E‘^ -F'^ £2 _ p2 

2Fp^^X^ 


£2 - £2 


( 2 . 21 ) 


Thus, we obtain the eigenvalue for down quark, 6 ^°"° by replacing p^ and a into 
their eigenvalues pz and u (= 0 , 1 , 2 ,---): 


4^" = \/(r ± \I-1QBvf + pI (for f = 1, 2,...) 

^7=0 = (for -r = 0) (2.22) 


with £3 = —£ and Q = —e/3. 

Thus, the single-particle energy of quarks with flavor f = u or d for up and 
down quarks can be expressed as 


e 


/ 

Pz,y,v 


/ 


F + r]^2\Qf\Bi 


+ Pz 


f z. = 0 , 1 , 2 ,--- 

\ i/ = l, 2,--- 


for 

for 


(2.23) 


with Qu = 2e/3, Qd = —e/3 and r/ = ±1. 

The thermodynamic potential 0 can be expressed in terms of the vacuum expec¬ 
tation values of the Hamiltonian H and the particle number operator N, in which 
the quarks occupy the energy levels from that with the lowest energy to that with 
the Fermi energy. Thus, we obtain ^ as 


= 3 ■ 1 E l(S"..r.+ - f) + (4,0,+ - M)1 

+3'7 E [(4,0.-- f) + (4,0,-- f)! + ^ . (2,24) 

Pz,!^ -</^) 


where V represents the volume under consideration and the factor 3 represents the 
color degree of freedom. Here, it should be noted that the sum with respect to pz 
and v can be regarded as the following R 1 : 


1 _ f dpz 

Pz,V 


\Qf\B 

2-k 


f(v) 

‘^max 



U=U 


f(v) 

min 


(2.25) 


Usually, i® replaced to f . Here, + Py can be regarded as 2\Q\Bv (= 

Px). Thus, the following correspondence may be understood: JJ dp^dpy — f 2'Kp±dp± = 
2^/ = 27r|Q|i3/di/ = 2n\Q\Bj:^. 
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Thus, the thermodynamic potential can be obtained finally as 



(2.26) 


where pp, and i^mai are determined by the condition ep^^u,r) < P- 
2.2. F = 0 case 

Next, let us consider the F = 0 case. For Q > 0 and Q < 0, Eos. (|2.lip and 
(j2.19h are valid with F = 0. Thus, the following equation for the two-component 
spinor cj) is obtained: 


-E + + \Q\B{2a}a -|- 1) 


^a^-\Q\B 


= 0 


(2.27) 


where the upper (lower) sign in front of corresponds to the case Q > 0 (Q < 0). 
Replacing p^ and a^a by their eigenvalues and cj) = *((/>i, (^> 2 ), the above equation is 
written as 


(—+ 2\Q\Bi'^ (j) I — 0 , 

2 

{-E^+pl + 2\Q\B{p + l))cj)^ =0. 

1 

Thus, the single-particle energy E is obtained as 

r ^pI + 2\Q\Bu 
\/ pz + 2\Q\B{iy + 1) 


(2.28) 


(2.29) 


This result is included in Ea. (j2.2.Sp . Thus, the thermodynamic potential Fq can be 
calculated as 


/ ' 

■ 


dp, \Qf\B 


r- „/ 


f=u,d 


27r 27r 


lv=0 


v=0 


+ ( a/pz + 2|Q/|7?(z^-I-1) — p 


dpz 

f=u,J-Po 


f=U,d'^~ 


\Qf\B{\Pz\ - p) 


dpz 

f=u4''-P° 


r _ 

2|Q/|R^ ^pl + 2\Qf\Bv-p 


y=l 


(2.30) 


where po and uL should be determined later. 
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§3. The thermodynamic potential in three cases: i^>/r, 0<T'</r and 

F = 0 


First, let us integrate out with respect to pz in the thermodynamic potential 
(j2.26p and (|2.30p . The condition (F + r/y^2|Qj|Fz^)^ + gives the range of 

integration with respect to pz- Namely, 

- + <Pz< (^F + rJ^J2\Qf\B^ (= pp) .(3.1) 


By using the following integration formula. 


ldpzV{F + X)^+pl = {pz + V{F + X)^+pI) , 

(3.2) 

we can carry out the integration with respect to pz in Ea. (l2.26l] . which leads to 


<p = — V 

27r ^ 


\Qf\B 

2tt 


Xv) 




( ! -\ ^ - {F + r]J2\Qf\Bvy 

+ {f + p^2\Qj\Bu\ In— V /V I /I ; 


F + pyj2\Qf\Bu 


+ 


2G 


f,M 


= S + A E mf\B Y. 9vimf\Bi^) , 


2G 87r2 


(3.3) 


/=w,d;77=± 


iv) 

U=Uf 

f,rn 


(3.4) 


{ti) iv) 

where Fj and Fj are determined by the condition 

-{f + r]^j2\Qf\Bi^ > 0 , 
which guarantees that pz is real. Here, Qrjix) is defined by 

drjix) = -p\lp? - (F + t/v^)2 + {F + pyGf In ^ ^ ■ (3.5) 

'' F + Py/X 

3.1. F > p case 

For r/ = 1, the condition (I3.4h is not satisfied for any F. On the other hand, 
for 7/ = — 1, if F < yj2\Qf\Bv, the condition (13.4p gives y/2\Qf\Bi' < ^ + F. If 
F > y/2\Qf\Biy, the condition (|3.4p gives y/2\Qf\Bu > F — p. Thus, we summarize 
the condition for v as follows: 




\{F-pr] 

+ 1 < p < 

liF + pf] 

[ mf\B \ 


[ 2\Qf\B \ 


- X Ti/T 


(3.6) 
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where [• • • ] represents the Gauss symbol. Thus, the thermodynamic potential can 
be expressed as 


_ V 

27r ^ 


f=u,d 


\Qf\B 

27r 


f,M 


E [-- (f'- \/2|<3 ,|Bi.) 

=^f7l 

+ (f - J2\Qf\Bu] In - \ ^ - - 

V V ; F-^2\Qf\Bv 


+ 


2G 


f,M 


Svr^ 


2\Qf\B Y. g.{2\Qj\Bv) + 


f=u,d 


u=u^-^ 

j,m 


2G 


(3.7) 


with ()3.6p . 

3.2. F < g, case 

For T] = 1, the condition (13.4p gives the condition 0 < ■\/2\QJ\Bi' < g, — F. 
On the other hand, for ry = — 1, if T < the condition (|3.4I1 gives F < 

< fJ, + F. If F > ^y2\Q flBu, the condition (|3.4h gives 0 < ^/2\Qf\Bu < 
F. Thus, we summarize the condition for ly as follows: 


I <i^< 


L 2\Qf\B 

if^ + Ff 


= .,(+) 


= V 




for ?y = 1 , 

for ?y = — 1 . 


2\Qf\B \ 

Thus, the thermodynamic potential is (|3.3p with (|3.8p for ry = ±. Namely, 


(3.8) 


(v) 

'^f,M 


■^‘“V^ + 5^2 E 2|<3/|b E 9,(2|'3/|b> 


(3,9) 


/=w,d;77=± 




3.3. F = 0 case 

Next, let us consider F = 0 case. Thus, the thermodynamic potential Fq can be 
given in ()2.3np . which we show again: 


Ft) = 




J 

M r 


+3 V r ^210,15V 


V=1 


\]pI + 2\Q f\Bv - g 
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4^2 E 


f=u,d 




E m\BY. 


f=u,d 


12=1 




Q Q IVl 

\Qf\B^,^ + J^ J2 2\Qf\BY,9omf\Bi.) , 

12=1 


471^ 


f=u,d 


f=u,d 


goix) = —g\/ g? — X + xln p— 


(3.10) 


with 


,/ - 


9 


[2\Qf\B\ • 


(3.11) 


§4. Approximation of the thermodynamic potential by replacing 
summation by integration with respect to the Landau level 


In the thermodynamic potential (12.2611 and (I2.30p . the quantum number v, which 
labels the Landau level, has to be summed up. However, since it is interesting 
to consider the spontaneous magnetization, it may be assumed that the external 
magnetic field B is small and finally B becomes 0. Therefore, let us replace the sum 
with respect to u by an integration approximately.^ 

In general, let us consider a function /(x). Here, we introduce a small quantity 
a and let us consider the Tailor expansion around x = av as follows: 


f a{u—l) 


dxf{x) = f 
J a 


dx 


a^u—l) 


f{au) + 


dx 


{x — au) + 


2 


(x — av)"^ + 


= af{aT^) - f'(au) + ^a^f''{au) + 

2 6 


(4.1) 


Thus, the following relations is obtained : 




E / 


rauM 

dxf{x^ dxf{x) 

J CLUrn 

=a ^ f{av) 

i2=!2m + l 


2 ’-'M 3 

Y + Y f'\o.v) + • • • . (4.2) 

I2=I2m + l ^2=^'„J + l 


Here, it should be noted that the dehnition of integral can be used when a is in¬ 
finitesimally small, namely, 

rauM 

a ^ f'{au) = / dxf{x) = f{auM) - f{ai^m) , (4.3) 

i2=i2,„+i 
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and so on. Thus, adding af{ai'm) on both sides of Ea. (j4.2ll . a useful approximate 
formula is obtained as follows: 


^ rai^M fj 

^ = / dxf{x) + - [/(oz^m) + fiavm)] - [f'{avM) - f'{avm)\ 

J aum. ^ ^ 


+ 


(4.4) 


After here, let us approximate the thermodynamic potential in the case of small B. 
For F > fi, Ea. (l3.7p is approximated by using the formula (14.41) as follows: 


^ ^ 3 

F = -1- 

2G 87r2 


f 

/ Pmax 

^ dpxp±g-{pi) 

J=U,a - ‘ min 

+ \Qf\B +5-(2|g/|i?i^^)) 

_2JQ^ (<7'-(2|Q;|Sz.|7^) - g'_{2\Qf\Bi.^-J^)) + • • • 


(4.5) 


with ([33]) . Here, = ^J2\Qf\Bu^ and pLax = ^J2\Qf\Bu^, respectively. 

For F < fi, the thermodynamic potential ()3.9p is approximated by using (14.21) 
directly as 


^ = S +A E 2iQ/i^ 


2G 87r2 


f=u,d 


t/(+) 

f,M 


f,M 


5 +( 0 ) + E g+{2\Qf\Bi^) + E 9-imf\Bi 


U=1 


U=1 


f=u,d 


3 r /"PmS . , 

Y1 ^ dpxp±g±{p'i) + \Qf\B (^g±{2\Qf\Biyylj) - 


f=u,d;r}=:t 


2\Qf\^B^ 


g'^{2\Qf\Bu^;^l)-g'^{0)) 


+ 


(4.6) 


where = [{p =F F)'^/{2\Qf\B)] and p^al = ■\j2\Qf\Bu^^lj, respectively. 
For F = 0, the thermodynamic potential (|3.1U|) is approximated as 


Fn = — 


dvr 2 


E i^/iv 


f=u,d 


3 r / 

+33 E / 9o{x)dx + \Qf\B \^go{afU-j^) - go 

X * . J d 0 


f=u,d 




[g’o{2\Qf\Bui,)-g'o{0))+--- 


(4.7) 
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Here, since B is small, u may be regarded as a continuum variables. Remem¬ 
bering the condition which determine s and namely, g±{2\Q= 0, 

and g±{0) / 0 with (|3.5j> . then Eas. (|4.5jl . (14.6p and (14.711 are simply calculated by 
performing the integration as 


IT'..3 rp2 


for F > fj. , 


F = ^ 


TT^ 


4 


F‘^g,\ Fg? 


arctan 


F\ g+yjg?- F2 

H-In-;;- 

12 F 


+ g + 0 (B^). 


v'm'^ - F 2 

for F < g , 


F = Fo = -^+0{B^) 


for F = 0 , 
where we used the integration formulae: 


(4.8) 


j dyyy/ g"^ - {F + gyf = g^ - {F + gyf 


-Y -^{F + gy) + ^{F + gyf 


Fg'^ 


arctan 


{F + gyfyjg^ - {F + gyf 
g^ - {F + gyY 


dyy{F + gyY hi 


- {F + gy^ 


F + gy 


= \/iF - [F + gyY 


-Y + ^{F + gy)- ^{F + gyf 


Fg^ 

H—— arctan 


{F + gy)^/g^ - {F + gyY 


-g^ + {F + gyY 


+ ^{F + gyf [-4F -h 3(F gy)] In 


h + VlF - {F + gyY 


F + gy 


(4.9) 


§5. Spontaneous magnetization for the quark matter under the 
tensor-type interaction between quarks 

First, let us calculate the spin polarization F with R = 0 as a function of the 
chemical potential g. This task has been done in our previous papers or IISI). 
The thermodynamic potential can be expressed as 

= j: Y_,.)+6.^ (F-M)+g. (5.1) 
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where the factor 6 represents the color and flavor degrees of freedom and pi = Px, 
P 2 = Py andp 3 = pz- Here, can be replaced to the integration f d^pj (27r)^. 

Integrating the three-momentum, we have obtained the thermodynamic potential as 


^B=o 




i?2 

1 

2G 

vr^ 

i?2 

Fp 

2G 

27r 


-^- {SF'^p, + 2p^) + Fp? arctan 

p + ^p?- F2 

-w in-- 

4 F 


F 


for F < p 

for F > p . 

(5.2) 


It should be here noted that Ea. (l5.2p is identical with (14.8p with B = 0. The gap 
equation is derived by d<pB=o/dF = 0. The quark number density pq is also derived 
through the thermodynamic relation pq = —dF/di^^ for the solution of the gap 
equation, F = Fyam- 


Pq 


= 


vr^ 


(^min + ) V P^ - + ^F^inP arctan 


7? 

A IT 


3G 


2 - ^2 . J 


p 


p 


for F < p 

for F > p . 

(5.3) 


Figured] shows (a) the spin polarization condensate F and (b) the baryon density of 
quark matter divided by the normal nuclear density, po = 0.17 fm“^, as a function 
of the quark chemical potential p. Here, we used a parameter G = 20 GeV“^ 



Fig. 1. (a) The spin polarization condensate F is depicted as a function of the quark chemical 

potential p. (b) The baryon number density of quark matter divided by the normal nuclear 
density po = 0.17 fm“® is shown as a function of the quark chemical potential. Here, pc ~ 0.407 
GeV is the value at which the spin polarization occurs. Also, in the region of /r > pp ~ 0.5605 
GeV, the spin polarization condensate is greater than the chemical potential, F > p. 
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which has been used E3 in our previous papers.®’® About = fie ^ 0.407 GeV 
which corresponds to the baryon density being 3.53 po, the spin polarization appears 
against the free quark phase. Of course, if the two-flavor color superconductivity is 
considered, it has been shown that the spin polarized phase appears about p = 0.442 
GeV {pB « 5.85po)-‘i^ 

We derive the spontaneous magnetization through the thermodynamic relation. 
The spontaneous magnetization per unit volume, A4, is defined by 


M 


d<P 


B=0 


(5.4) 


In the right-hand side, 5 = 0 is adopted. However, from (|4.8I) . the thermodynamic 
potential does not depend on B linearly. Thus, the spontaneous magnetization is 
equal to zero under this consideration, namely 


M 


d<P 


= 0 , 

B=0 


(5.5) 


even if the quark-spin polarization occurs, while the magnetic polarization which is 
proportional to B may appear. 


§6. Spontaneous magnetization originated from the anomalous 
magnetic moments of quarks 

In the previous section it has been found that no spontaneous magnetization 
of polarized high density quark matter is predicted by the normal coupling to an 
external magnetic field. However, the spin polarization F should be observed in 
some way. We know that the quark has an anomalous magnetic moment. In this 
section, let us investigate the effect of the anomalous magnetic moment of quarks. 

The effect of the anomalous magnetic moment is introduced at the level of 
the mean field approximation in this paper. Here, it is known that the anomalous 
magnetic moments reveal ad® 

^ ) , pu = 1.85 pN , Pd = -0.97 pN , 

Mrf J 

PN = = 3.15 X 10"^^ GeV/T . (6.1) 

2mp 

We introduce the effects of the anomalous magnetic moment in the Lagrangian den¬ 
sity within the mean field approximation as 

Ca = C- , (6.2) 

If the vacuum polarization is taken into account, the parameter G with rather small value 
gives the same results quantitatively such as G = 11.1 GeV“^ under the standard three-momentum 
cutoff yl = 0.631 GeV, although we do not consider the vacuum polarization in this paper. 
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where = dfj,Ay — dyA^ and F 12 = = —B. Here, L is nothing but Ea. (|2.3p . 

We only take /r and i/as/r=l, v = 2 and /U = 2, v = 1 because the magnetic field 
has only 2 ;-component. Then, in the mean field approximation, Ea. (l2.3|) is recast 
into 


Ca = iipYDfj,ip - i^iFsTs + ijaB1)F3'iI; . (6.3) 

Here, F 3 T 3 = FI where 1 is the 2x2 identity matrix for the isospin space, which are 
denoted in Eas. (l2.2h and (I2.9p . Thus, we introduce the flavor-dependent variables 
Ff as 

Ff = F + fifB, namely Fu = F + ^uB , Fd = F + ^aB . (6.4) 

Therefore, the Lagrangian density can be expressed as 

Ca = - Tp{F + ijlaB)F: 3 iI) - — 

F^ 

= - V’FT'sV' - ^ . (6.5) 

Thus, we learn that the variable F should be replaced lo Fj = F + njB for each 
flavor except for the last term in (16.5p . In this replacement, we can derive the 
thermodynamic potential in the same way developed in §4. Namely, the results 
including the effects of the anomalous magnetic moment of quarks should be obtained 
by replacing F into Fj for each flavor in the previous calculation developed in §4, 
except for the term originated from the mean field approximation, F‘^f2G. The 
results are summarized as follows: 


^ = -*-> = 4 S ^ + S + tor r>A.. 

f=u,d 



Thus, the spontaneous magnetization per unit volume, A4, can be derived thor¬ 
ough Ea. ()5.4p . that is 


M 


d<P 


B=0 


E 

f=u,d 


Wf 


B=0 


dJl 

dB ' 


(6.7) 


From (16.6p , we can derive the spontaneous magnetization originated form the anoma¬ 
lous magnetic moment and the spin polarization: 


M 




— (/^ti + l^d) 1 


B=0 


(6.8a) 
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U [GeV] 


Fig. 2. The spontaneous magnetization per unit volume M is depicted as a function of the quark 
chemical potential /r. The vertical axis represents the magnetization Af in a logarithmic scale. 
For < 0.407 GeV, the magnetization does not occur because the spin polarization does not 
appear, F = 0. 


M = -Y. 


80 , 


f=u,d 

1 

27r2 


dFf 


B=0 


dB 


2Ffi -\/arctan 
X il^u F Fd) 


F 




-F^ In 




F 


~ 2G 


(6.8b) 


where the gap equation d0/dF = Ylf=ud^^/^^f + d{F‘^/{2G))/dF = 0 is used 
from the second line to the third line in ()6.8bl) . It should be noted that if = 0, 
the spontaneous magnetization disappears because F = 0 in Ea. ()6.8bp . 

In FiglJl the spontaneous magnetization is depicted as a function of the quark 
chemical potential. The spontaneous magnetization is shown in SI unit. For n < 
fic ~ 0.407 GeV, the magnetization does not occur because the spin polarization does 
not appear, F = 0. At /r = /ic, the spontaneous magnetization suddenly appears. 

§7. Magnetic field of hybrid compact star 

As seen in the previous section, for ^ > Hc, the spontaneous magnetization in 
the quark matter appears. Here, the spontaneous magnetization per unit volume 
Ai may be regarded as the magnetic dipole moment. Under this identification, we 
can calculate the strength of the magnetic filed yielded by the spontaneous magne¬ 
tization. As is well known in the classical electromagnetism, the magnetic field at 
the position r, namely magnetic flux density B{r), created by the magnetic dipole 
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Magnetc flux density 

Fig. 3. The magnetic flux density is shown in a logarithmic scale as a function of the quark chemical 
potential /i in the case (a) = 1 km where the quark matter occupies 0.1 % of the total volume 

of star ((47rrq/3)/(47ri?^/3) = 0.001) and (b) = 2.15 km where the quark matter occupies 1 

% of the total volume of star ((47rrq/3)/(47ri?®/3) = 0.01), respectively. It should be noted that 
the SI unit, Tesla, is converted to Gauss, namely, 1 T (Tesla) = 10“* G (Gauss). 


moment m can be expressed as 


dvr 



3r{m ■ r) 

1^5 


(7.1) 


where /to represents the vacuum permeability. In our case, m = (0,0, Ad x y), where 
V represents a volume because Ad is nothing but the magnetization per unit volume. 

Here, let us consider the hybrid star with quark matter in the core of neutron 
star. Let us assume that the hybrid (neutron) star has radius i? = 10 km. If there 
exists quark matter in the inner core of star from the center to Vq km, the strength 
of the magnetic flux density on the surface at the north or south pole of hybrid star 
is roughly estimated as 


B, 


Ho f Ad ‘iz^M\ 
dvr V R^ ) 



2Mrl 


[T] . 


(7.2) 


Figure[3](a) and (b) shows the magnetic flux density as a function of the quark chem¬ 
ical potential h the case (a) = 1 km and (b) Vq = 2.15 km where quark matter 

occupies (a) 0.1 % and (b) 1 % of the total volume of star ((d7rrq/3)/(d7ri?^/3) = 
(a) 0.001 and (b) 0.01), respectively. It should be noted that the SI unit, Tesla, is 
converted to Gauss, namely, 1 T (Tesla) = lO'^ G (Gauss). In the magnetar, the 
strength of the magnetic field at the surface of star is near 10^^ G. Thus, in our cal¬ 
culation, if quark matter exists and the spin polarization occurs, a strong magnetic 
field about 10^^ or 10^^ Gauss may be created. 


§8. Summary and concluding remarks 


,/i[GeV] 


It has been shown that spontaneous magnetization occurs due to the anomalous 
magnetic moments of quarks in the high-density quark matter under the tensor-type 
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four-point interaction. In the Nambu-Jona-Lasinio model as an effective model of 
QCD, the tensor-type four-point interaction has been introduced. As for the tensor- 
type four-point interaction, this interaction term was also introduced to investigate 
the meson spectroscopy, especially, for vector and axial-vector mesons.^^ Owing to 
this interaction, the spin polarized condensate appears for each quark flavor in the 
region of large quark chemical potential. It has been shown that the spin polarized 
condensate leads to spontaneous magnetization of quark matter due to the anoma¬ 
lous magnetic moments of quarks. Also, it has been pointed out that spontaneous 
magnetization does not occur if there exists no anomalous magnetic moments of 
quarks. 

In this paper, furthermore, the implication to the strong magnetic field in com¬ 
pact stars such as the hybrid stars has been discussed. If there exists quark matter 
in the core of neutron stars and then the quark number density is rather high, spon¬ 
taneous magnetization may occurs. If quark matter occupies a volume with 1% of 
neutron star, the strength of magnetic field at the surface of neutron star is of order 
of 10^^ Gauss, which is comparable to the strength of magnetic field in the so-called 
magnetar. 

As indicated in this paper, the spin polarized condensate appears in the region 
with a large quark chemical potential. Thus, the chiral symmetry is broken in this 
model because the starting Lagrangian density is constructed with chiral symmetry. 
Under the magnetic field, the chiral symmetry is broken and it is shown that the 
chiral condensate grows at most linearly as a function of the magnetic held 
Further, in Ref J23p . the tensor-type four-point interaction was introduced in the NJL 
model in the case of one quark-havor. In that paper, both the chiral condensate and 
spin polarized condensate are equally treated at hnite temperature system. Thus, it 
may be an interesting problem that the coexistence of the chiral condensate and the 
spin polarized condensate is also considered at hnite baryon density system. 

As for the implication to compact stars, it has been shown that there is a possi¬ 
bility of the existence of the massive hybrid quark stars with two solar mass under 
the strong magnetic held.® Thus, the investigation of the equation of state of quark 
matter in the spin polarized phase revealing the spontaneous magnetization may be 
one of the interesting future problems. 
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